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The entanglement in a pure state of TV qudits (d- dimensional distinguishable quantum particles) 
can be characterized by specifying how entangled its subsystems are. A generally mixed subsystem 
of m qudits is obtained by tracing over the other N — m qudits. We examine the entanglement 
in the space of mixed states of m qudits. We show that for a typical pure state of N qudits, its 
subsystems smaller than N/3 qudits will have a positive partial transpose and hence are separable or 
bound entangled. Additionally, our numerical results show that the probability of finding entangled 
subsystems smaller than N/3 falls exponentially in the dimension of the Hilbert space. The bulk 
of pure state Hilbert space thus consists of highly entangled states with multipartite entanglement 
encompassing at least a third of the qudits in the pure state. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Lx 



I. INTRODUCTION 



and also extend the analysis to any dimension of qudits. 



Quantum information is a rapidly developing field ex- 
ploiting the peculiar entanglement properties of quantum 
states. Its applications can be broadly divided into two 
general types: quantum computation {|, ||, and quantum 
communication Quantum communication (including 
quantum cryptography and quantum teleportation) can 
be framed in terms of repeated use of pairs of entan- 
gled qubits - the entanglement properties of two qubits 
have been well characterized, and a number of analytical 
measures of entanglement are known j|, |[ ||, 0, [|, |). 
However, for multi-qubit systems, which are essential for 
quantum computation, few entanglement measures can 
be calculated even for pure states. Despite these difficul- 
ties, arrays of qubits have been the focus of recent atten- 
tion |o[ fyj |§, ||, H, H), though often only pairwise 
entanglement has been considered in such systems. The 
more general question of entanglement in multi-qudit sys- 
tems (d-dimcnsional quantum particles) is also impor- 
tant since most real systems (e. g., atoms) have more 
than two states. 

In this paper we tackle the question of entanglement 
in larger qudit systems by deriving a bound on the size 
of typical entangled mixed states that are subsystems of 
larger pure states. Previous numerical work by some of us 
p6[ suggested that qubit subsystems somewhat less than 
half the size of the pure state, typically have no usable 
entanglement. In other words, since pure states are typi- 
cally highly entangled, the entanglement is distributed in 
multi-partite entanglement involving around half or more 
of the qubits. This work was limited by computational 
resources to small numbers of qubits (less than 14) . Here 
we provide analytical bounds for any number of qubits, 



We emphasize that we are interested in average prop- 
erties of all possible pure states. It is easy to construct 
states with entanglement properties that lie outside our 
bounds. For example, the so-called W states (symmetric 
superposition of one qubit in state |1) with the rest in 
state |0)) have only pairwise entanglement |l5| for any 
value of N. However, such states are of small measure 
and do not contribute significantly to the average prop- 
erties of pure states sampled from the whole of Hilbert 
space. 



To further motivate our investigations, we note a con- 
nection to the the question of the origin of the speedup 
in quantum computing. Recently Jozsa and Linden in 
|]l7f argue that for an exponential speedup over classi- 
cal computation, the state of a quantum computer using 
pure states will have multi-partite entanglement encom- 
passing an arbitrarily large fraction of the total number of 
qubits in the computation. However, they also point out 
that the entanglement cannot be inferred to be the cause 
of the speedup, rather, it is one of a number of properties 
of the typical states in the Hilbert space of the quantum 
computer that are necessary for such a speed up. Our 
work confirms this property of typical states in Hilbert 
space, and also provides an analytical lower bound on the 
number of parties in the multi-partite entanglement for 
a given size of Hilbert space. 



The paper is organized as follows. First, in Sec. [n], we 
introduce clear definitions of the system we are consider- 
ing and the entanglement measures we will use to analyze 
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it. Next in Sec. [II, we derive the analytic bounds, fol 
lowed by comparison with numerical results in Sec. 
We finish with conclusions in Sec. [v|. 



IV 
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II. DEFINITIONS AND ENTANGLEMENT 
MEASURES 

In a finite dimensional Hilbert space Ti there exists a 
natural, unitarily invariant, Fubini-Study measure ptFS, 
induced by the Haar measure on the unitary group. A 
system of N qudits has a Hilbert space of size d N . We 
will investigate entanglement in such systems by sam- 
pling randomly from the set of all pure iV-qudit states 
according to the natural measure (ips and calculating 
their average, which is to say, typical properties. This will 
give us useful information about the majority of states in 
Hilbert space. Having chosen a random pure state of N 
qudits, we can divide it into two subsystems by partial 
tracing. A subsystem of size to is obtained by partial 
tracing over the other N — m qudits. This partitioning 
could represent the system of interest plus the environ- 
ment, for example, or two different parts of one system. 
From the mathematical point of view one obtains in this 
way a certain probability measure pN,m in the space of 
mixed quantum states acting in a d m dimensional Hilbert 
space induced by the natural measure pps on the space 
of pure states in Ti d N [jl8| . This is a useful way to sample 
randomly from mixed states, since unlike pfs for pure 
states, there is no unique natural measure for the space 
of mixed states. 



A. Entropy of subsystems 

Consider now a pure state quantum system described 
by a density operator p = \<p)(4>\ acting in a composite 
Hilbert space Hmk = 'Hm^'Hk of dimension MK. The 
reduced density matrices defined by the partial trace, 
Pm '■= trKP and pk '■= tiMP characterize both sub- 
systems. It is well-known that the von Neuman en- 
tropies of both subsystems are equal, Sm — S(pm) — 
— trpMhipA/ = S(pk) — Sk- The value of Sm shows 
how entangled the two subsystems are with each other. 
If Sm = 0, there is no entanglement and the composite 
state may be factorized, p — pm ® Pk, but this is not a 
typical case. The mean entropy of a subsystem averaged 
over the natural measure /x^s is given by 
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where M < K and the approximation holds for <C 
M < K. This result was first conjectured by Page fl9[| 
and later proved in j2(], 21. 

In our case the subsystems consist of to and N — to 
qudits, so substituting M = d m and K — d N ~ m into 
(jj) we obtain the mean entropy (S)N,m of the subsystem 
analyzed, where we are now using the subscripts to re- 
mind us of the number of qudits in the pure state (N) as 
well as the number in the subsystem (to). Equation (|l|) 
tells us that on average, the smaller subsystem has nearly 
maximal entropy, showing that the two subsystems are 



highly entangled with each other, a typical pure state is 
highly entangled, see also |||. The particular value of 
Sn.tu provides a good measure of entanglement between 
the two subsystems, but says little about entanglement 
between the qudits within a single subsystem. Asymptot- 
ically, for (JV — > 00), to < N/2 will tend to be maximally 
mixed and therefore separable, but to obtain useful re- 
sults for finite N, we need further information on the 
size of the region around the maximally mixed state that 
contains separable states. 



B. Purity and mixedness 

Next we will characterize how mixed a typical m-qudit 
state p m might be. A convenient measure is the purity 
p3| defined by 



r[p m ) := tr(p^). 



(2) 



Pure states are defined to have r(p m ) = 1 while mixed 
states have l/d m < r(p m ) < 1. We will also refer to 
R := l/tr(/4) = l/R as the inverse participation ratio 
(IPR) of the mixed state p m . The IPR thus ranges from 
1 (pure) to d m (maximally mixed) , and larger IPR means 
the states are more mixed, while larger purity means the 
states are more nearly pure. 

Consider random states drawn according to the natu- 
ral measure pfs on the space of pure states in a MK 
dimensional Hilbert space. As before, we apply a partial 
trace to obtain the reduced density matrices pm and px- 
In refs. [^8[ |i| it was shown that the average purity 
of Pm and pk is equal to 



M + K 
MK + 1 



(3) 



This result does not depend on the particular way in 
which the Mif-dimensional Hilbert space is decomposed, 
but only on the initial and the final dimensionality of 
the spaces. Thus it also holds in the problem we are 
analyzing here of pure states of ./V-qudits reduced to m— 
qudits by partial tracing. Substituting AI — d m and 
K = d N ~ m the average purity of the system reads 



jN-r 



(r) 



N.m 



d N + l 



(4) 



In other words, the averaging has been performed with 
respect to the induced measure pN,m- 

Numerical investigation of the induced measures shows 
that for to < N/2 the probability distributions P(r) 
are concentrated close to the mean value (r). This is 
illustrated in fig. [j], where P(r) is plotted for d = m = 2 
and N = 3 to 6 in the main plot, and for d = 2, to = 3, 
N = 5 to 10 in the inset. The scale of the purity axis 
on the inset has been adjusted (the scale is in fact linear 
in R = 1/r) to show the data near the minimum purity 
more clearly. 



3 




purity, r 

FIG. 1: The numerically obtained distribution of the purity 
r for (main graph) mixed states of m = 2 qubits obtained by 
partial tracing of random pure states of (right to left) N = 3, 
4, 5 and 6 qubits, and (inset) mixed states of m — 3 qubits 
for N = 5 to 10. The vertical dashed lines indicate the values 
of rppT from eq. (ph. 



it has been shown || that all the states sufficiently close 
to the maximally mixed state, with R > d 2 — 1 have a 
positive partial transpose (PPT). This condition defines 
the maximal ball inscribed in the convex body of mixed 
states. Very recently it was shown that all mixed states of 
a bipartite system belonging to the maximal ball posses 
not only the PPT property but are also separable p0[ . 
The former result may be generalized for a multipartite 
case: any state p m (within the d m dimensional Hilbert 
space) for which 



R > Rppt = d m —l so that r < r PPT 



1 



(5) 



has the PPT property with respect to any possible trans- 
positions of the subsystems. We give the proof of the 
generalization in Appendix [A|. 



III. BOUNDS ON ENTANGLEMENT IN 
SUBSYSTEMS 



Fixing the number to, of the qudits in the final system 
and varying the initial number N of qudits involved in 
the pure state we get a sequence of probability measures 
which sweep the space of the mixed states close to the 
manifold of a constant purity. Hence it is justified to use 
the mean value (r) as a good estimate of the actual value 
r for a typical subsystem. 



C. Peres criterion, positive partial transpose 

We will now describe a simple test for entanglement in 
multiparticle mixed states due to Peres p6|j . Consider a 
multiparticle mixed state p acting in a d m dimensional 
Hilbert space, which describes the system consisting of to 
qudits. An 771-particle system can be split into two parts 
containing j and m — j particles respectively. The partial 
transpose is applied to the j particle part of the system 
producing p Tj := (1l ® Tj)p. Here Tj denotes the trans- 
pose operation in the d 3 dimensional subspace, while the 
dimension of the subspace that has not been transposed 
is L = d m ~i . We say that the TO-qudit state has the 
PPT property, if p Tj is positive for all possible values 
of j (it is sufficient to take 1 < j < to/2) and choices 
of j particles, i. e., with respect to all of its possible 
splitting into two subsystems. Conversely, a state which 
fails this test is described as being NPT (negative partial 
transpose). NPT states have some useful entanglement 
that could be distilled with some probability into one 
or more maximally entangled states of two qubits. This 
test does not distinguish bound entangled states |27|, |||] 
from separable states, but since bound entanglement is 
relatively rare |29|], and for most purposes it is the free 
entanglement that is useful, this still provides a useful 
characterization of entanglement. 

For bipartite systems with Hilbert space dimension d 2 , 



We are now ready to provide an analytical relationship 
between the size of a subsystem and its typical entangle- 
ment properties. 



A. Bound on PPT region 

Using the fact that the distributions of the purity P(r) 
are narrow for the regions of interest (to < N/2), and 
thus (r) jv,m is a good estimate of the actual value of r in 
a typical subsystem, we apply this to the critical value 
rppT obtained in the previous section. If the mean value 
(r) averaged over a certain induced measure is smaller 
that the critical value rppx, (or (R) > Rppt), the ma- 
jority of random states distributed with respect to this 
measure are localized inside the maximal ball and are 
PPT. We therefore combine eqs. (||) and (|) to calculate 
a relationship between N and to such that we expect 
almost all subsystems of to qudits will have a positive 
partial transpose. We find (r) N,m < Tpt when 



d N > d im {\ - dr m - d-' Im ). 



(6) 



Taking the logarithm (base d) gives us an estimate of how 
many qudits the initial pure state should contain, such 
that after the reduction to m-qudits the obtained mixed 
state is PPT on average, 



A^ppt > 3to + log d (l - dr m - dr 2m ). 



(7) 



The log d term is always negative, becoming rapidly 
smaller for either m > 2 or d > 2. Hence for systems 
initially consisting of 



N > 3to 



(8) 



qudits we expect the subsystems of size m to be PPT 
mixed states with a considerable probability. For the 
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simplest case of m = d = 2, the smallest integer num- 
ber larger than the right hand side of (Q) gives Ap PT = 
int[6 — log 2 (16/ll)] + l = 6. For a fixed number N this es- 
tablishes a bound on to below which subsystems of this 
size have on average no useful entanglement in an N— 
qudit pure state. This is the main result in this paper. 

It is also possible to calculate the maximum value of 
the von Neuman entropy for states on the surface of the 
max imal ball, see Appendix [b|. This value S c , given in eq. 
(Bl), can then be combined with eq. (]]]) for the average 
entropy of a subsystem, to obtain a similar relationship 
between N and to. The result is (as it should be) the 
same as eq. (^) , but with the restriction to ^> from the 
approximate form of eq. ([!]) . Our derivation as presented 
using eqs. (||) and (^) is applicable to all < to < N. 



B. Estimation of the transition region 

The above result is a bound on the average properties 
of typical states that identifies the range of parameter 
values (N > 3to) where there is a high probability that a 
subsystem of size m is PPT. We would also like to charac- 
terize the range of N, m for which with high probability 
the subsystems of size to are entangled. This region is 
not simply the inverse of the PPT region because while 
the criterion we are using in eq. (Q) is sufficient for ensur- 
ing the subsystem is PPT, it is not a necessary condition. 
There are many states which are PPT that have r > rppx 
right up to separable pure states with r = 1. 

Since we do not, in general, know the size of the region 
of PPT states outside the maximal ball defined by eq. 
(Q), our approach is to identify a region of the body of 
mixed states of a considerable size, which contains en- 
tangled states only. There is a one parameter family of 
mixed states (known as the generalized Werner states) 
containing to qudits defined as a mixture of the maxi- 
mally entangled pure state Y$>) and the maximally mixed 
state Ijm j 



(l-e) 



I d m+e\*)(V\ 



where 



(9) 



(10) 



and e is a real parameter between and 1 specifying 
the proportions of the mixture. These states have been 
shown to be entangled for |BL, H 



£ £ p/n.t 
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(11) 



and strictly separable for e < e ent (there are no bound 
entangled states in this family). This boundary at e ent 
defines a set of entangled states of positive measure ^] 
as can be seen as follows. The states p e lie on a line that 



is on an axis of rotational symmetry in the set of mixed 
states. This line intersects the set of separable states at 
£ — Sent ■ The set of separable states is convex, so beyond 
the hyperplane normal to the p e line there exists a set of 
entangled states of a positive measure (see also 34 ). 

It is straight forward to show that the purity for p e at 
this boundary point (e = £ en t) is given by 



2d 



^ent 



(d m + dy 



(12) 



Thus there exists a greater-than-zero probability to en- 
counter entangled states with r > r en t. However, we 
don't know that entangled states actually dominate over 
PPT states for this value of r. Since Werner states are 
among the most entangled states J3E 36, 57) for a given 



purity r, we might actually guess that this estimate will 
turn out to be too tight. Nonetheless, it provides a useful 
independent check on our previously derived bound, so 
we present it anyway. 

Equating (r) M,m from eq. (^) with r ent from JT^ ) al- 
lows us to establish an estimate for the entangled side of 
the transition between entangled and PPT states. It is 
easily shown that (r)jvm > r ent when 



Nent < 3TO-2+log rf [l + (2d+l)d- m +(d+2)d 



l-2ml 



(13) 



For all d and to, the log d term is positive and tending to 
zero for increasing d and to. Hence for all finite N with 



N < 3m — 2 



(14) 



we expect the probability of finding the subsystem of 
size to entangled (not PPT), to be non-zero. Looking 
at the smallest case [Q, d = to = 2, we have N ent = 
int[4 + log 2 (ll/4)] = 5. Let us emphasize again, that 
this is an estimation only; for this (or lower) values of 
N the probability of finding entangled subsystems of size 
to is non-zero, but it does not rule out the existence of 
substantial numbers of PPT subsystems of size to, so the 
transition might occur effectively for even smaller values 
of N. 

Now let us fix the initial size of the pure states of N 
qudits while we decrease the size of the final system of to 
qudits. For to = N the probability of finding a separable 
pure state is equal to zero. As m is reduced, the states 
obtained by partial tracing over N — to qudits become 
increasingly mixed, and the probability of finding PPT 
states increases. Putting both results together we may 
characterize quantitatively a transition region 



3m - 2 < N < 3r 



(15) 



which does not depend on the qudit size d, in which we 
estimate that PPT subsystems come to dominate over 
entangled subsystems. (Note that the inequalities have 
reversed because we are identifying the region where the 
ratio of PPT to NPT (negative partial transpose) subsys- 
tems of size to is of order unity, rather than the regions 
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FIG. 2: Numerical results showing the transition from entan- 
gled to PPT subsystems. For each m, points are shown for 
which Pppt G (1%,99%) for qubits (■) and qutrits (grey v)- 
The solid/dashed lines represent the range in eq. (D-M. 
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FIG. 3: Probability of finding that a subsystem of m qudits is 
entangled in random pure states sampled uniformly over the 
Haar measure for qubits for m = 2 (■) and m = 3 (filled o). 



where one case dominates over the other.) These are 
surprisingly tight estimates of the boundary between en- 
tangled and PPT states in terms of the number of qudits 
in the state. However, it must be remembered that the 
number of qudits is a logarithmic measure of the size of 
the Hilbert space. The addition of one qudit enlarges the 
Hilbert space by a factor of d. 



IV. NUMERICAL COMPARISONS 

In |l6[ ], for the qubit case d = 2, it was shown numer- 
ically that entanglement in a mixed m-qubit subsystem 
of a typical iV-qubit pure state falls off sharply towards 
zero with increasing N for values of to < 5 and N < 13. 
These numerical results suggest the transition from NPT 
states dominating to PPT dominanting is complete by 
N ~ 2to + 3, which falls outside the range ofeq. fllfj) 
for to > 4. This is at the edge of the results presented in 
ijbSf , so we ran further numerical studies to confirm this 
divergence between numerical and analytical results and 
also, since our analytical results apply for any dimension 
of the constituent quantum particles d, to test qutrits 
and confirm they follow the same pattern. 

Using straightforward C programs optimized for effi- 
ciency to generate the large data sets required for useful 
statistics, we generated random pure states in the d N 
dimensional Hilbert space, traced over N — to qudits, 
and analyzed whether the remaining mixed state, with a 
<i m -dimensional Hilbert space, had a positive or negative 
partial transpose (PPT or NPT). In fig. || we represent 
the transition region in the N-m plane for qubits up to 
to = 6, N — 15 and qutrits up to to = 3, N — 8, combin- 
ing data from flfj] with our new data for N — 14 to 15 
(qubits) covering the transition for to = 6, and qutrits 



for N — 4 to 8. Our new results confirm the trend. For 
to > 4, the transition from NPT to PPT for a given sub- 
system size to is completed for N < 3m — 2. This is due 
to our estimates not taking into account the full size of 
the set of PPT states, being based only on the maximal 
ball inscribed in this set. 

There is one more feature of this transition that re- 
quires comment. As noted in the transition is very 
sharp, the proportion of PPT subsystems of a given size 
to becomes essentially zero to numerical accuracy with 
the addition of just a few qudits to the pure state. It is 
not obvious a priori that the transition should be sharp. 
One might reasonably expect it to be smoother, with a 
typical pure state of N qudits containing finite amounts 
of to— partite entanglement for all 2 < to < N. In order 
to obtain such a sharp transition, the proportion of NPT 
subsystems of given size to must fall exponentially in the 
size of the Hilbert space d N of the pure state. This is hard 
to check numerically, but the data we have do support 
this, see fig. |[ This sharp transition and exponential fall 
off allows us to assert that this bound on the multi-partite 
nature of entanglement in typical states of Hilbert space 
corresponds to the distinction made in |l7j] between sets 
of states that can be simulated efficiently classically, and 
sets of states that cannot. It also justifies our statement 
in the introduction that exceptions to our bound (^) do 
not contribute significantly to the average properties of 
typical pure states. 



V. DISCUSSION AND CONCLUSIONS 

We have obtained analytically a lower bound on the 
size of subsystem of a typical pure state that can be 
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expected to contain any entanglement within the sub- 
system. For a typical pure state of N qudits, its sub- 
systems of size m < N/3 can be expected to be PPT, 
i.e. to have no usable entanglement within them. We 
then checked numerically and confirmed that the tran- 
sition from NPT to PPT subsystems actually occurs at 
significantly larger m nearer to N/2, which is the asymp- 
totic limit as N — > oo. Our results emphasize the im- 
portance of multipartite entanglement in quantum in- 
formation processing. The types of states necessary for 
exploring the main bulk of Hilbcrt space with quantum 
computers are highly entangled, but the entanglement 
will not be evident if only a small subset of the qudits 
are examined. For a quantum computer operating in a 
pure state, at least a third of the qubits would need to 
be processed to detect the entanglement present. 



operations Tj of partial transpose, so the state p is PPT. 
□ . 

It is illuminating to disc uss a simple geometric inter- 
pretation of the condition (|Al[) , which defines a certain 
subset of the convex body of mixed states M. acting in 
n = d m dimensional Hilbert space. The center of the 
n 2 — 1 dimensional set A4 is given by the maximally mixed 
state p* = I/n. Introducing the Hilbert-Schmidt metric 



D 



(pi,p 2 ) = (tr( Pl - P2 ) 2 ) 1/2 (A2) 



HS 



it is not difficult to compute the distance of any state p 
with eigenvalues {xj} to the center of M., 
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D Hs(PiP*) = J2{ a 
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APPENDIX A: PPT PROPERTY FOR MIXED 
STATES OF MULTI-PARTITE SYSTEMS 

In this appendix we prove the following 
Proposition. Any mixed state p of m qudits which 
satisfies the condition 



R(p) 



1 

trp 2 



1 



(Al) 



has the PPT property, i. e., all its partially transposed 

T- 

matrices p„l are positive. 
We start the proof by invoking an algebraic 
Lemma. Let A be anx n non-zero Hermitian matrix 

and consider the real number a := tiA/ (tiA 2 ) 1 / 2 . If 

a > yjn-l then A > 0; 

the proof of which is given by Mehta, |58| , inequality 

9.21, p. 217. Let us apply it for the analyzed density 

matrix partially transposed in an arbitrary way, A = 

p Tj , which remains a Hermitian matrix of size n = d m . 

Any operation of partial transpose Tj does not change 



the traces, so trpU = trp = 1 and tr(p 



1/R. Therefore, the coefficient a 2 — 1/trp 2 is just equal 
to the inverse particip atio n ratio R. If a mixed state p 
satisfies the condition (Al), the assumption a > \Jn — 1 
is fulfilled. Hence matrices p Tj are positive for all possible 



Hence the condition R =const determines a set of points 
equidistant from p». 

The boundary of M. is defined by the condition detp = 
0, so any state belonging to it has at least one eigen- 
value equal to zero. The state p a with the spectrum 
•f— ^r, — . . . , —^r, 01 is closest to p* among the states 

L n — l 5 n— 1 1 ' n — 1 ' J o 

belonging to the boundary of M. Its inverse partici- 
pation ratio, R(p a ) = n — l, coincides with the critical 



value from the constraint (Al). Thus all states satisfying 
it form the maximal ball inscribed into the convex body 
of mixed states M. and centered at p*. 



APPENDIX B: VON NEUMANN ENTROPY ON 
THE BOUNDARY OF THE MAXIMAL BALL 

In this appendix we outline the proof of the following 
Proposition. If the von Neuman entropy of a state 
p is larger than S c , then p belongs to the maximal ball 
and its inverse participation ratio R is larger than n—l 
where n is the dimension of the Hilbert space of p and S c 
is given by 



S c = In n 



-ln2+^ ^ln 

n n 



1 



n-2 



(Bl) 



In the asymptotic limit n — > oo, Sc — * ran, i.e. maximal. 
State p a , with the spectrum 



Xn 



1 



1 



n—l n—l 



.0 , 



(B2) 



belongs to the boundary of the set of mixed states and 
to the maximal ball inscribed in it since 



R = 



1 



tr(p 2 ) 



= n — 1, 



(B3) 



see Appendix Its entropy S(p a ) — ln(n — 1) is the 
smallest among all states located at the surface of the 
hyperspehere R = n — 1. To show this it suffices to per- 
turb the spectrum keeping R constant and compute the 



7 



entropy variation. Due to the symmetry of the problem, 
the state pi, which belongs to this hypersphere and has 
the largest entropy, is located opposite to p a with respect 
to the maximally mixed state p* = I/n and we have 



(Pb + Po)/2 = P* 



This gives the spectrum of pb 
n-2 



n-2 



n(n — 1) ' ' n(n — 1) 



(B4) 



(B5) 



which evidently belongs to the maximal ball R - 



Computing its entropy S(pb) 
value S c given by (Bl). □. 



one obtains the critical 
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